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A Test of The Source Galerkin Method 

D. Petrov a , P. Emirdag a , G. S. Guralnik a 

a Physics Department, Brown University, 
Providence, PJ 02912 

Some results of the ongoing development of our Source Galerkin (SG) nonperturbative approach to numerically 
solving Quantum Field theories are presented. This technique has the potential to be much faster than Monte 
Carlo methods. SG uses known symmetries and theoretical properties of a theory. In order to test this approach, 
we applied it to (f> 4 theory in zero dimensions. This model has been extensively studied and has a known set 
of exact solutions. This allows us to broaden the understanding of various properties of the SG method and to 
develop techniques necessary for the successful application of this method to more sophisticated theories. 



1. Introduction. 

The Source Galerkin method is being devel- 
oped as a flexible alternative to Monte Carlo ap- 
proaches to solving Quantum Field Theories. 

1.1. Overview of the Method 

To illustrate the application of the Source 
Galerkin approach to solving Quantum 

Field Theories we consider a theory defined by 
a Lagrangian C. All information of the theory is 
contained in the generating functional 



Z(j) = J P0exp J {-C+ 3 4>)d d x 



(1) 



This satisfies a Schwinger-Dyson set of differential 
equations. A procedure to determine an approx- 
imate solution of these equations is specified by 
our SG technique. To start, a trial solution or 
an ansatz must be introduced. In general, it is 
constructed from a predefined set of trial func- 
tions. The proper choice of trial functions allows 
us to take advantage of symmetries of the theory 
as well as of other known analytical properties of 
the generating functional. A simple example of 
an approximate solution is: 



Z„ 



exp 



jxG X yjy~\~ I jwjxHwxyzjyjz ~f" 



•(2) 



A set of residues is obtained by substituting this 
ansatz into the Schwinger-Dyson equations. The 
undefined coefficients of the trial solution are de- 
termined by solving a system of nonlinear equa- 



tions obtained by projecting the residues on the 
trial functions and requiring these projections to 
be equal to zero. The procedure outlined above 
guarantees that the error associated with the ap- 
proximate solution converges to zero in the mean 
as the number of members in the set of trial func- 
tions goes to infinity. 

1.2. Motivation 

In order to demonstrate the validity of the 
method, it was applied to the 0(3) nonlinear er 
model Q. This model is asymptotically free and 
is a useful toy model for approaching Non-Abelian 
gauge theories. First order solutions were ob- 
tained, but during our attempts to extend this 
work to higher orders it became evident that a 
better understanding of the general properties of 
the SG method is necessary. Theorems guarantee 
that Galerkin methods produce an approximate 
solution, which converges to the exact solution as 
the number of terms goes to infinity. In practice, 
the approximate solution can only have a fairly 
small number of parameters. Therefore, we must 
ensure that high accuracy and rapid convergence 
can reasonably be achieved with only a few terms 
in the ansatz. The fact that the error approaches 
zero in the mean implies that the accuracy of the 
result is heavily affected by the choice of inner 
product. We attempt to investigate the effective- 
ness of the Source Galerkin approach by using it 
to solve </> 4 theory in zero dimensions. We study 
its performance for several choices of trial func- 
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tions and scalar products. 

2. Ultralocal Model 

After introducing an external source j, we de- 
fine the ultralocal model by the following la- 
grangian 



(3) 



Theoretical solutions for this theory can be ob- 
tained by expressing generating functional as a 
power series 



k=0 

The coefficients of this series are given by O 
.U(n,t) + (-l) n pU(n,-t) 



(4) 



a-2n 



(2n - !)!!■ 



2n». 

«2n+l = r\~ t) 



U(0,t)+pU(0,-t) 
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V(±,i) U(0,t) + pU(0,~t)- 



Here U and V are parabolic cylinder functions, 
coefficients p and a fix boundary conditions. 

2.1. Numerical Solution 

In analogy with the theoretical solution, a trun- 
cated power series of order N can be chosen as 
a trial solution. After obtaining the residues by 
substituting this expression into the Schwinger- 
Dyson equations, we need to eliminate source de- 
pendence by projecting them on the trial func- 
tions. Two possible ways to define scalar product 
are presented below. 



j n R(j)dj, n 0...Y 



and 



PR(j) exp 



dj, n = 0..N 



(5) 



(6) 



Both definitions are localized in the region close 
to j = since eventually we want to set the source 
to zero in order to compute physical values. Us- 
ing these definitions produces similar results as 
shown in figure 1. However, the second expres- 
sion is more general and it is easier to extend 
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Figure 1. Order by order relative error is plot- 
ted for N=5 and residual equations derived from 
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J j n R{j)dj with respect to the range of integra- 
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to calculations with space-time dimensions. In- 
spection of these results shows that the error in 
determination of coefficients dj increases rapidly 
at high orders. In spite of this, the generating 
functional can be determined with accuracy as 
high as several parts in 10 7 for a significant range 
of j. However, it must be noted that correct de- 
termination of the proper range of integration or 
a value of the parameter e is crucial for achiev- 
ing high accuracy. The SG approach by itself 
does not provide an algorithm for setting these 
parameters. In a real problem when the exact re- 
sult is not available it is necessary to have some 
external source of information which allows us to 
determine the correct values of the parameters. 



2.2. Solution in Terms of Hermite Polyno- 
mials 

A set of Hermite polynomials can be defined by 



.2 i cTexp(-^} 
W-.li) = i-l)"«>x 1 .^| ^ -■ (7) 
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Figure 2. Order by order relative error is plotted 
for N — 7 and N = 5. Residual equations derived 



from / H n (j)R(j)exp 

— oo 

the parameter e. 
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dj with respect to 



These polynomials are orthonormal under the fol- 
lowing inner product 



H k (x)H l (x)exp 



dx = S k 



(8) 



This suggests a way to improve our implemen- 
tation of the Source Galerkin procedure for the 
ultralocal theory. Trial functions can be con- 
structed as a linear combination of Hermite poly- 
nomials of different orders, while the expression 
shown above is used to define the scalar prod- 
uct. This choice leads to significant simplification 
of the system of equations which give the coeffi- 
cients in the ansatz. Figure 2 shows the depen- 
dence of relative error in determination of values 
at on the parameter e. From this graph we ob- 
serve that there is a significant range of values of 
the parameters in which numerical results remain 
stable. We can use this fact to resolve the prob- 
lem outlined in the end of previous section. In a 
real problem, it would be possible to set all the 
parameters without any knowledge of theoretical 
solution just by finding solutions for all possible 
values of parameters and then choosing a solution 



which is stable within the required accuracy for 
some range of the parameters. The accuracy of 
determination of the generating functional using 
this method is consistent with the power series 
ansatz approach. This implies that it is limited 
only by error introduced by the method used to 
solve the system of equations produced by the 
Source Galerkin procedure and by the numerical 
precision of the hardware used for computation. 

3. Conclusion 

This investigation demonstrates that very high 
precision computations can be performed using 
the Source Galerkin procedure. On the other 
hand, the method is very sensitive to the choice 
of parameters. Application of this method to a 
more general set of problems will require devel- 
opment of some consistent method of determina- 
tion of proper values of parameters. Current re- 
sults suggest that use of orthogonal trial functions 
may simplify this task. Source Galerkin technique 
is still under development. Additional work is 
needed to fully investigate it's properties and to 
apply it to various field theories. 

4. Acknowledgements: 

We wish to thank R. Easther, D. Ferrante, 
Z. Guralnik and S. Hahn for helpful discus- 
sions. This research was supported in part by 
DOE grant DOE grant DE-FG02/19ER40688- 
(Task D). 

REFERENCES 

1. P. Emirdag, R. Easther, G.S. Guralnik, S.C. 
Hahn, New Numerical Methods for Quantum 
Field Theories on the Continuum. Nucl. Phys. 
Proc. Suppl. 83 (2000) 938-940 

2. S. Garcia, Z. Guralnik, G.S. Guralnik, Theta 
Vacua and Boundary Conditions of the 



Schwinger-Dyson Equations, bcp-th/9612079 



